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THE SMITH CLOUD: HIGH- VELOCITY ACCRETION AND DARK-MATTER CONFINEMENT 
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ABSTRACT 

The Smith Cloud is a massive system of metal-poor neutral and ionized gas (M gas > 2 x 10 6 M ) 
that is presently moving at high velocity (Vqsr ~ 300 km s _1 ) with respect to the Galaxy at a 
distance of 12 kpc from the Sun. The kinematics of the cloud's co metary tail indicates t hat the gas 
is in the process of accretion onto the Galaxy, as first discussed by ILockman et a l. ( 2008). Here, we 
re-investigate the cloud's orbit by considering the possibility that the cloud is confined by a dark 
matter halo. This is required for the cloud to survive its passage through the Galactic corona. We 
consider three possible models for the dark matter halo (NFW, Einasto, Burkert) including the effects 
of tidal disruption and ram-pressure stripping during the cloud's infall onto and passage through the 
Galactic disk. For the NFW and Einasto Q dark-matter models, we are able to determine reasonable 
initial conditions for the Smith Cloud, although this is only marginally possible with the Burkert 
model. For all three models, the progenitor had an initial (gas +dark matter) m ass that was an order 
of magnitude higher than inferred today. In agreement with ILockman et al.l . the cloud appears to 
have punched through the disk rs 70 Myr ago. For our most successful models, the baryon to dark 
matter ratio is fairly constant during an orbital period but drops by a factor of 2 — 5 after transiting 
the disk. The cloud appears to have only marginally survived its transit, and is unlikely to retain its 
integrity during the next transit w 30 Myr from now. 
Subject headings: 



1. INTRODUCTION 

The Smith Cloud is perhaps the best known and most 
studied member of the high- velocity clouds (HVCs). 
While covering over a third of the sky, little is known 
about the nature of these clouds, with many suggestions 
as to their origin. While we now know the distances to 
a number of HVCs, their orbits are m ore difficult to de- 
termine. Since its discovery in 119631 . the Smith Cloud 
was originally believe d to be an extension of the Galac- 
tic disk (|Smitb]!l963ft but it was later shown to be sep- 
arate d kinematically from the Galactic disk (|Lockmanl 
I1984D leading to its subsequent classification as an HVC 
(jWakker fe van Woerdel3ll997ft . 

While the origin of the cloud is still unknown, with 
suggestions of gas belonging t o an infalling dwarf galaxy 
(Bland-Ha wthorn et alJ 11998ft or gas lifted from the 
Galactic disk (jSofue et al.ll2004l ). the distance, mass and 
velocity of the cloud are now well established. Dis- 
tance indicators based on stellar absorption-line brack- 
eting, Ha flux, and the kinematic distance of gas inter- 
acting with th e cloud, provide a consistent estimate of 
12.4± 1.3 kpc (ILockman et al.ll200llPutman et alJl2003t 
iWakker et al.ll2008l). The HI mass of the cloud is found to 
be 1 x 10 6 M (jLockman et al.ll2008ft . with WHAM data 
indicating a comparable or greater amo unt in ionized 
hydrogen (jHill. Haffner fc Reynolds! 12009ft. The velocity 
comp onents have also been determined (|Lockman et all 
2008) allowing the orbit of the cloud to be calculated 



for the first time. The Smith Cloud's importance is that 
it provides insight into gas at a more advanced stage of 
accretion than the Magellanic Stream. 

Earlier studies have shown how easi ly gas disrupts 
as it passes through a hot corona flMoore fc ~ 
1994tlQuilis fc Moorell2q01l : lBland-Hawthorn et all 



Heitsch fc Putm an 2009) . It is therefore remarkable that 
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1 The Ei nasto profile is sometimes referre d to as the Sersic pro- 
file (e.g. Zait, Hoffman, & Shlosman 2008). These have the same 
functional form but since the Sersic profile arises from a projected 
distribution, we prefer to use a radial function. We note, however, 
that the projection of the E inasto function is nontrivial analytically 
UMazure. fc C apclato 2002). 



a massive, high-velocity gaseous system has survived so 
close to the disk. In light of recent developments, we 
consider the prospect of a confining dark halo surround- 
in g at least so me HVCs to be plausible, as first proposed 
bv lBlitz et all (|2001ft . and therefore worthy of further in- 
vestigation. A number of ultra-faint dwa rf galaxies have 
now been discovered in t he Galactic halo ( Simon fc Gehal 
2007; Ki rbv et al J 120081 IGeha et al.ll2009ft . What is par- 
ticularly striking is that the inferred dark halo mass ex- 
ceeds ~ 10 7 Mq ev en for systems with very few stars 
(Stri gari et all 12008ft . These observations hint at the 
prospect of dark matter dominated galaxies containing 
HI that have not experienced star formation. Complex 
H is one such 'dark galaxy' candidate, with a mass domi- 
nated by dark matter if it is self-gravitating (|Simon et al.l 
l20M ). 

Here we suggest that the Smith Cloud is another candi- 
date for a dark galaxy, illustrated in Fig. [TJ with the gas 
partly stabilised by a dark matter component. We calcu- 
late the subhalo required for the survival of a gas cloud 
with properties similar to those observed in the Smith 
Cloud today. In $2]we discuss the properties of three pos- 
sible dark matter profiles that may encompass the Smith 
Cloud. In f}3] we discuss the orbit of the Smith Cloud 
within a realistic Galactic potential. In SJH we present 
the equations that are used to calculate the ram pres- 
sure stripping of the subhalos along their orbit. In §5\ 
we present the results of the semi-analytic models and 
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Fig, 1. — Hi image of the Smith Cloud (Fig. 1. in Lockm an et alj 
(2008)) with superimposed core and halo for three different dark 
matter models (NFW, Einasto, Burkert) discussed further in J5] 
The tidally stripped mass today is 1 — 2 X 10 8 Mq in dark matter. 
The core size assumes an initial gas mass at apogalacticon of 1.1 X 
10 7 Mq, although most of this gas is distributed along the orbit 
of the Smith Cloud today. We note that the Burkert profile has 
no surviving core, and the NFW core is displayed in white for 
visibility. 

discuss the implications of these results for the Smith 
Cloud. In SjS]we summarize these results and discuss the 
limitations of our model. 

2. DARK MATTER PROFILES 

Numerical simulations of dark matter halos indicate 
that they tend to be characterized by two parameters, 
typically a radial scale length and a density. One of the 
first of these dark matter profil es was inspired by rotation 
curves in dwarf spiral galaxies (|Burkertlll995f ). the 'Burk- 
ert' profile, which incorporates a core of uniform density 
followed by a steady decline with increasing radius. Early 
dissipationless N-body simulations however revealed that 
the dark matter core is more cusp-like. The most widely 
us ed of the cuspy profiles is the 'N FW profile introduced 
bv lNavarro. Frenk fc White! (|1997f ) . But more recently, it 
seems that neither of these profiles are sufficiently flex- 
ible to describe what is now being seen in simulations 
involving 10 10 particles, with a more appropriate mo del 
provided by the Einasto profile (Springc l et al.ll2008l ). 

All of these profiles may be described by a density func- 
tion fp(x) with the density at any point given by p(x) = 
Pdsfp(x), where p ds is the characteristic density. Simi- 
larly the mass may be described by M(x) = Md s f m (x) 
and the potential by ip(x,v s ) — v^f ip (x), where v s is the 
characteristic velocity of the halo. The gas density for a 
hydrostatic cloud is given by f gas (x, v s ,c g ) = nn / n>H 



gas* 

c is 



,0 



the isothermal sound 



exp[-(v^/c 2 g )f v (x)], where w, 
speed of the gas, and nu,o is the nuclear gas density. 

These profiles are defined in terms of a scaled radius 
x = r/r s , wher e r„ is a scale radius. These param eters 
are defined by ^Sternberg. McKee fc Wolfirell2002f ): 



M ds = -np ds r s 
pds = Ap„ 
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Properties of DM subhalo (M vir =1.2x10 a M ) 
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Fig. 2. — Comparison of the different subhalo density and mass 
profiles considered by this work prior to tidal stripping (r B = 
1.04 kpc). 
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where A is the overdensity of the dark matter ~ 340 at 
z = 0, p u = 2.8 x 10 -30 g cm -3 is the average density 
of the universe, M V i r is the virial mass, and x V i r is the 
scaled virial radius. x V i r is correlated, with a standard 
deviation a, to the virial mass in the virial mass range of 
10 s - 10 11 M by (ISternberg. McKee fc Wolfir^[2002ll : 
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Convenie nt forms for the NFW and Burkert profiles 
are given in ISternberg. McKee & Wolfirel (puul ). In Ta- 
ble 5 of that paper, we note that there is a typograph- 
ical error in the NFW f gas equation which should read 
f gas = e -3{vs/c g ) 2 ^ + x j3(v„/c g ) 2 /x_ In Table [U we give 
our formalism for the Einasto profile, and include the 
Burkert and (corrected) NFW profiles for completeness. 
The three density and mass profiles are shown in Fig. [2] 
for median subhalos of virial mass 1.2 x 10 9 M , with 
corresponding scale radius r s = 1.04 kpc and profile de- 
pendent characteristics density determined by $2$. 

For dark matter subhalos moving within the Galac- 
tic environment, the radial profile will be altered in two 
competing ways. The condensation of baryons contracts 
the inner core of the subhalo while ram-pressure strip- 
ping leads to an expansion of the subhalo. The overall 
impact of ad iabatic contraction can be calculated itera- 
tively using (jBlumenthal et al.lll986f) : 



r [M b (r) + M DM {r)] - r l M l (r l ), 



(G) 



where M (r) is the final baryon mass distribution, 
M^uif) is the final dark matter distribution and Mi(ri) 
is the initial total mass distribution at an initial orbital 
radius r^. By assuming the orbits of baryonic material 
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f p (x) = x-' L (l + x)^ 
/„,(*)= 3 + 
/».(*,«.)= 3 [l-^=£±=l] 


/ p (x) = (l + a; )- 1 (l + ^ 2 )- 1 

/«.(*) " I [ ln ^+ x2) + ln(l + x) tan" 1 s] 

= §[(! + |) tan" 1 x - (1 + |) ln(l + x) 

+ l(l-I)ln(l + ^)] 
/ gos (a:,« s ,c s ) = [ e -(l+V*)ten l I(1 + J ,)(H-l/«) 

(1 + x 2)(l/2)(l/x-l)](3/2)(« s /c g ) 2 


/ p = exp[-2/«(x« - l)]/4 

/m = /3 7 (3/a,2x a /a) 

f v (x,Vs) = /3[2 1 ^ a a- 1 ^ a y(2/a,2x a /a) 

-j(3/a,2x a /a)/x - 1] 
fgas(x,v s ,c g ) = exp(-v%/c*f lp (x)) 

P = (3/4)8- 1 /«e 2 / Q, a- 1+3 / a 



TABLE 1 

NFW, Burkert and Einasto profiles, where 7 is the lower inc omplete gamma function, x = r / r s is t he scale radius, v s is 

THE HALO CIRCULAR VELOCITY AND Cg IS THE GAS SOUND SPEED; CF. STERNBERG. McKEE & WOLFIRE (2002), TABLE 5. THE FOUR 
QUANTITIES IN EACH COLUMN ARE THE DENSITY PROFILE f p , THE DARK MATTER MASS PROFILE / m , THE DARK MATTER POTENTIAL PROFILE 

ftp, AND THE GAS DENSITY PROFILE fgas- 



do not cross, then Mdm{t) = (1 — F)Mj(rj) where F is 
the baryon to dark matter fraction. We have investigated 
the impact of adiabatic contraction and find that it has 
a negligible influence in modelling the Smith Cloud. 

We note that the virialization time of the dark matter 
halo can also be important. Here the overdensity of the 
dark matter compared t o the average density is given by 
(jBrvan k Normanlll998T ): 



18tt 2 + 82a; - 39x 2 

~ n(z) ' 

where x = Q(z) — 1, and for a flat universe 

fi m (l + z) 3 



(7) 



(8) 



The physical density of a halo is directly proportional to 
this overdensity factor such that pas = PoA/A where p a 
is the central density of a halo that virialized at z = 0. 
This factor also contributes to other halo properties such 
as the scale radius r s oc A -1 / 3 and the scale velocity 
v s oc A 1 / 6 . 

3. MODEL SETUP 

We consider two models of evolution, one in which the 
Smith Cloud is infalling for the first time, hereafter the 
Infalling Orbit Models, and a second model where the 
Smith Cloud has already been maximally stripped due 
to previous orbits, hereafter the Repeated Orbit Mod- 
els. These both share common features: (i) they have 
the same trajectory today, (ii) the dark matter halo has 
been tidally stripped down from some larger initial mass 
(M v ii) in an identical fashion before our calculations 
commence at apogalacticon. The important distinction 
is tidal stripping of the gas is possible in the Infalling 
Orbit Models but not in the Repeated Orbit Models; in 
both cases, ram pressure stripping by the hot halo is im- 
portant. For each case, the evolution of the Smith Cloud 
is considered for the NFW, Einasto and Burkert models. 

The evolution of the model clouds was calculated as 
a function of three variables: the initial virial mass at 
the time of formation (i.e. before the dark-matter halo 
fell into the Galaxy), the dark matter profile at this 
time, and the initial hydrogen gas mass at apogalacti- 
con. For both the Repeated Orbit Model and Infalling 
Orbit Models, the evolution of 7503 model clouds were 
calculated, corresponding to 61 logarithmically spaced 
virial masses in the range M v ; r = 5x 10 7 — 5 x 10 10 M© 



-150 Myr 




Fig. 3. — T he orbit of the Smi th Cloud, calculated using the 
potential from Wolfirc ct al. (1995). The current position is repre- 
sented by an unfilled circle and the Smith Cloud is travelling in the 
direction of the arrows, with heights below the disk represented by 
a dotted line. The Sun's position is shown as a filled circle on the 
Solar Circle. The thin dotted line represents the projection of the 
Smith Cloud onto the disk. The disk is represented by a solid line 
at 30 kpc. 



and 41 logarithmically spaced gas masses in the range 
6 - 1 x 10 8 Mr. 



M gas = 1 x 10 



i . 



The orbit of the Smith Cloud was calculated us i ng the 
position and velocity data from lLockman et alj (2008) 
for the tip of the Smith Cloud: (R, z) = (7.6, -2.9) kpc, 
(vr,v^,v z ) = (94,270,73) k m s" 1 . The form o f the 
Galactic potential is given by IWolfire et all (1995) nor- 
malized by a circular velocity of v c = 220 km s _1 at 
the Solar Circle. In Fig. [31 we show the predicted orbit 
of the cloud system. In agreement with Lockman ct al 
(2008), we find that the Smith Cloud has intersected the 
disk ~ 70 Myr ago and will pass through the disk again 
in ~ 30 Myr. 

For all subhalo models, we investigate the effects of dy- 
namical friction on t he orbit trajectory. Th e formalism 
used is described by Uiang fc Binnevl (|2000l ): we point 
out that the value for the circular velocity in their Ta- 
ble 1 should be v c — 235 km s _1 (not 181 km s _1 as 
quoted) to be consistent with their analysis. But over 
the past few hundred million years, dynamical friction is 
found to have only minimal effect, even in the high mass 
limit. This is because, once again, the impact of gas 
loss from the subhalo close to the disk is found to dom- 
inate the evolution of the subhalo. We assume that any 
drag between the model clouds and the Galactic corona 
is negligible and does not affect the orbit. 

Each model cloud is considered to be a dark matter 
potential well filled with gas in isothermal hydrostatic 
equilibrium. We assume a primordial helium abundance 
nHe/nn = 1/12 and metallicity of Z/Zq = 0.1. We 
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also assume the gas has a temperature of 1.2 x 10 4 K 
and adopt an ionization fraction of 50% for the Smith 
Cloud , slightly below the new l y upd ated H + /H° ra- 
tio in Hill . Haffner fc Reynold! (|2009| ). This tempera- 
ture and ionization fraction then give a sound speed of 
c g = 11 km s . The gas is distributed in the potential 
well according to the gas density profile riH{x,v s ,c g ) = 
nH,ofgas(x,v s ,c g ), where x = r/r s is the scale ra- 
dius, v s is the halo circular velocity given by Q and 
fgas{x, v s ,c g ) is given in Table [TJ 

For the initial dimension of the model clouds, the sound 
crossing time is 200 Myr, falling to about 30 Myr at the 
disk. We therefore begin each orbit at the apogalac- 
ticon, approximately 17 kpc from the Galactic centre, 
which took place 150 Myr ago, and f» 100 Myr after 
the last transit of the disk. We hence assume the cloud 
has regained hydrostatic equilibrium by this point and is 
the basis for our initial assumption of hydrostatic equi- 
librium. Due to the comparatively large sound crossing 
time we assume that the radial gas density profile does 
not change throughout the calculation, only being trun- 
cated, however, due to the lower sound crossing time at 
the disk we use a more restrictive condition as the cloud 
transits the disk (Q. 

In order to trace the long-term impact of ram-pressure 
stripping, we must define the density profile of the hot 
coronal halo and the Galactic gas disk. The accelera- 
tion due to the stellar density profile is already included 
in the Wolfire potential. The Galactic corona is as- 
sumed to be a mix of H+He at Th = 2 x 10 6 K with 
th e density calculated from the gravitational potential 
of I Wolfire et alj (|1995[ ). normalized to a plasma density 
of n ffn = 10~ 4 cm" 3 at 55 kpc from the Galactic Cen- 
tre (|Bland-Hawthornl 120081 ) . The density profile of the 
HI disk was calculated using an exponential drop-off in 
density both in the radi al and vertical dire c tion. We use 
the new prescription in iKalberla &: Pedes! ((2008) which 
includes the vertical flaring of the outer disk. 



4. ANALYTIC PROCEDURE FOR CLOUD/HALO 
EVOLUTION 

The Smith Cloud was evolved through a procedure 
that iteratively calculates the competing processes on 
the dark matter confined cloud. The initial radius of the 
gas comes from integratin g the gas density profile over 
space giving (modified from Sternberg, McKee & Wolfire 
[20021) : 



(|McCarthv et al.ll2008h : 



M gas (0) = A-KrlmHfgasixo)- 



P(0) 



583fc h T, 



fc J c Jo 



X fgas(.X^)dx : 

( 9 ) 

where r s is the scale radius of the halo, Xq = r(0)/r s , 
fgas(x) = f g as{x, v s ,c g ) is the gas distribution profile for 
constant scale velocity and sound speed. Here we have 
replaced Uh,o with Pf (nk(,T)f gas (xo) by assuming that 
initially the outer edge of the halo is bound by thermal 
pressure from the Galactic halo. 

At each step we calculate whether the gas at the cloud 
edge is prone to stripping by examining if the force 
applied by the external medium exceeds the maximum 
gravitational force in the subhalo at a given radius. For 
an isothermal spherical cloud, gas will be stripped if 



, nGM(R)p c 



2R 



(10) 



with p ga i and p c being the density of the Galactic gas 
and the outer subhalo gas respectively, v c the relative 
velocity between these two gases, and M(R) the amount 
of mass within a radius R. The mass within a radius 
R comprises the dark matter subhalo Mom and the gas 
within that radius M gas . 

If this stripping criterion is satisfied, then a 
shock passes through the gas, at a speed given by 
(|McCarthv et all 12008) : 



V fs 



Pgal 
Pc 



(11) 



This gas assumed to be no longer bound and is removed 
from the cloud. By assuming that the shock removes 
a shell of gas, the outer radius of the gas changes at a 
rate given by r — —Vf s /2. As Vf s tends to be much larger 
than the gas sound speed, c g , we assume that the density 
profile of the gas does not change, becoming truncated 
at the edge of the shock. 

For Infalling Models, the tidal stripping must also be 
considered. At each step, the tidal radius and tidal mass 
of the su bhalo are also calcu lated. The tidal radius is 
given by (jHavashi et al.ll2003f ): 



M(r t ) 



2 - 



R DM, 



gal 



M gal (R) OR 



MggljR) 

R 3 ' 



(12) 



where R is the distance between the centre of the sub- 
halo and the Galactic centre. Any gas outside of the 
tidal radius r t will also be removed, with the new gas ra- 
dius being the minimum of the tidal radius or the radius 
resulting from ram pressure stripping. 

The total mass can then be calculated by modifying 
© such that 

P(0) f x ^ 

M gas (t) = 4nr s m H f gas (x ) i 5g3fc T j x 2 f gas (x)dx, 

(13) 

where x(t) — r(t)/r s . The central column density is 
calculated using 



N„Jt) 



P(0) 



1.583fc h T r 



fgas(x)dx. (14) 



We iterate through the orbit in 1 Myr steps to determine 
the evolution of radius, mass and column density for the 
cloud/halo system. Convergence was tested by iterating 
through on 0.5 Myr steps for some surviving clouds, with 
the smaller step size producing less than 1% difference 
in mass calculations throughout the orbit. 

A special step, however, is considered when the Smith 
Cloud is moving through the Galactic disk. Consider a 
cloud with a gas density p c , radius r c moving at speed 
v c through the Galactic halo. Consider a gaseous disk 
where, at some radius, it has a local gas density pd and 
a half-scale height hd- Will the cloud survive its transit 
through the disk? 

We write down two shock conditions: one for the shock 
driven into the cloud by the disk, and one for the shock 
driven by the cloud into the disk. The cloud survives if 
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the shock through the cloud takes longer to propagate 
than the shock through the disk. The shock propagates 
into the disk gas at speed 



4:V r 



V s ,d 



3(1 + y/pd/Pc) 
and another into the cloud, at speed 



3(1 + y/p^/p~ c ) 



(15) 



(16) 



The shock needs to be fast enough to heat the disk gas 
significantly. 

The key to punching a hole is that the shock must go 
right through the disk gas before passing through the 
cloud. In other words we need roughly r c /v StC > hd/v Si d 
or equivalently v St d/v StC > hd/r c . This leads to 



V Pc/Pd > h d /r c 



(17) 



It follows that, to first order, the cloud will survive if 



N H ,c > —N Htd , 



(18) 



where Nh,c and Nhm are the column densities of 
the cloud and the disk along the orbit respectively. 
This is the kinetic e nergy argument presented by 
iBland-Hawthornl (|2008). We note that the dark mat- 
ter subhalo itself has negligible impact on the Galactic 
gas. 

We therefore expect that all gas in the Smith Cloud 
at a column density less than hd/r c ■ Nn,d is stripped 
when the cloud moves through the disk. As the passage 
through the disk occurs on a short timescale < 1 Myr, 
we treat this removal of gas as instantaneous, and in ad- 
dition to the normal ram pressure stripping that occurs 
during that time step. We also set our model disk to ro- 
tate with equal angular velocity to the Smith Cloud as it 
passes through at w 13 kpc, and a static disk otherwise. 
Although this choice is unphysical and in both cases con- 
tradict, the observed relative velocity between the Smith 
Cloud and the ISM, the cloud spends only a short time 
near the disk, and hence any additional gas stripped via 
ram pressure will contribute only a small correction to 
the mass lost from the cloud near the disk. 

5. MODEL RESULTS 

In our simulations, a large number of model clouds sur- 
vive to the present day, although not all of these are sim- 
ilar to the Smith Cloud as observed today. We interpret 
our results graphically through contour plots of proper- 
ties derived from the mass, radius and column density 
of the model clouds. The simplest properties — mass, 
column density — are displayed for the Repeated Orbit 
Model NFW profile in Fig. 0J By removing model re- 
sults which are greatly dissimilar to those observed for 
the Smith Cloud, we can find an allowed parameter space 
for the initial cloud properties (Fig. [5|). In Figs. 4-7, we 
present our results in terms of the mass of the tidally- 
stripped subhalo rather than the mass of the initial in- 
falling halo. 

The lower limit of the allowed parameter space is found 
by requiring the amount of gas in the model clouds to 
match the observations. The Smith Cloud has an HI 



mass of at least 10 6 M Q (iLqckman et all I2008D, and 
at least as much again in HII (Hill, Haff ner fc Reynolds! 
12009ft . giving a final gas mass of at least 2 x 10 6 M©. The 
excluded region corresponding to lower cloud gas masses 
is shown in green (Fig. [5]) bounded by a solid black line. 

We include several upper bounds to the allowed pa- 
rameter space. The universal baryonic to dark matter 
mass fraction is f2b/Sl m w 0.17 from the latest microwave 
background observations. 1 While there may be galaxies 
that exceed this limit (e.g. dwarf galaxies that form in 
tidal tails of strong mergers) , it is reasonable to assume 
that on average most dark matter subhalos do not ac- 
crete more than this fraction. Thus, we can probably 
exclude all clouds with baryon to dark matter mass frac- 
tions above this limit, displayed by the red shaded region 
bounded by a dot-dash line in Fig. El 

In dwarf galaxies, the highest HI column densities 
seen outside of star-forming regions are always less than 
2 x 10 21 cm -2 . This critical column density is found in 
regions with < la of optical emission (B. Ekta 2009, pri- 
vate communication). Since no stars are observed in the 
Smith Cloud, we include this upper limit to the HI col- 
umn density, shown by the light blue region bounded by 
a dotted line in Fig. EJ 

An extensive literature demonstrates that star for- 
mation is intimate ly associated with the pre sence of 
molecular gas (e.g. iCalzetti fc Kennicuttl 120091 ). For a 
canonical HVC metallicity of about 0.1 Z , it is un- 
likely that molecular gas can be measured directly from 
the Smith Cloud wi th existing submillimetre telescopes 
(|Ritcher et al.ll200"U ). This sets a stronger, but presum- 
ably less rigorous constraint on the allowed regions of the 
Smith Cloud; the excluded regions are shown in purple 
bounded by a dot-dash line in Fig. [5] The detection 
of molecular gas in the Smith Cloud would alter this 
limit from an upper limit, to a lower limit and hence 
alter the allowed parameter space to exist only in the 
purple region. All of the upper limit bounds shown in 
Fig. [5] have a strong dependence on the gas distribution 
within the model. Here we have assumed that the gas 
was initially in hydrostatic equilibrium with the profile 
only being truncated and not altered further; shocks that 
substantially modify the profile are associated with gas 
that gets stripped. 

The NFW and Einasto profiles restrict the allowed re- 
gions of the cloud to a narrow range of initial gas masses 
and subhalo masses, with the subhalo mass range re- 
stricted to 10 8 — 10 9 M Q and gas mass initially 5 x 10 6 — 
5 x 10 7 Mq, although we note that the allowed gas mass 
range is much larger at smaller subhalo masses than at 
the higher mass end. By comparison to the gas mass 
present today, the Smith Cloud will have lost between 
50% and 95% of the gas present at apogalacticon. A 
large fraction of this gas loss in all model clouds is due 
to the transit of the disk suggesting, suggesting that the 
Smith Cloud may not survive the next passage through 
the disk. The tidal stripping of dark matter does not sub- 
stantially affect the gas contained in the model clouds, 
with the gas contained within the protected core of the 
potential well. In the Infalling Orbit models, the Smith 
Cloud will also have lost approximately 50% of the dark 

1 Throughout the paper, we assume a ACDM cosmology with 
parameters (a s ,h, Si b , fi m , fi A ) = (0.8, 0.7, 0.04, 0.3, 0.7). 
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matter mass as well to tidal stripping, although most of 
this mass loss will have occurred before the semi-analytic 
calculations were begun. This scenario of an HVC in- 
falling for the first time, possibly after interacting with 
a third body, such as a dwarf galaxy, has the potential 
to increase the baryonic to dark matter ratio above the 
universal fraction, particularly if the orbit avoided dense 
sections of the disk. The Burkert profile retains a lower 
column of gas, due to a lower central density, with the 
effect that the region is not as well constrained as in the 
case of a NFW or a Einasto profile. It is this lower col- 
umn that leads to a higher mass requirement in order 
to survive traversal of the disk, which is dictated by the 
column density of the cloud. 

We display the mass history of three Infalling Orbit 
Model clouds corresponding to the same input parame- 
ters for each dark-matter profile (M v i r = 1.2 x 10 9 M 
[M tidal « 3x 10 s M ], M gaa = l.l x 10* Mq) in Fig. E The 
hydrostatic pressure bound gas radius is located within 
the tidally protected core of the dark matter subhalo, and 
therefore the mass history of the Repeated Orbit Model 
with the same initial conditions will be similar. 

As the subhalo loses dark matter to tidal stripping 
and gas to ram pressure stripping at different rates, the 
baryon to dark matter fraction will alter over the course 
of the orbit. This change, displayed in Fig. [7] for the 
Infalling Orbit Model clouds corresponding to the input 
parameters (M vir = 1.2 x 10 9 M [M tid ai ~ 3 x 10 s M ], 
M g as=l-1 x 10 7 M©) is particularly apparent at the tran- 
sit of the disk where a large amount of gas is lost but 
the dark matter component remains unchanged. Noting 
that the dark matter halo will have already been tidally 
stripped at the beginning of the evolutionary model sug- 
gests that HVCs may undergo radical changes in their 
baryon to dark matter fraction along their orbits. 

6. SUMMARY 

We have modelled the Smith Cloud by assuming that 
it is supported against total disruption by a dark-matter 



subhalo. By considering the evolution of this system over 
the past 150 Myr, we are able to compare our predicted 
properties to the known properties of the Smith Cloud to- 
day. Under the assumption of dark-matter confinement, 
the subhalo encompassing the Smith Cloud has a tidal 
mass of ~ 3 x 10 8 M at the present time. Even before 
this subhalo arrived at its present orbit, it will have lost 
50 — 90% of its dark matter and a comparable amount of 
progenitor gas due to tidal and ram pressure stripping by 
the Galaxy. In this respect, the likelihood of an infalling 
(dark-matter confined) "Smith Cloud" event is compa- 
rable to the rate of dwarf galaxy infall, i.e. roughly one 
event per Gyr. 

If the NFW and Einasto profiles are appropriate de- 
scriptions of dark matter halos, we find that the Smith 
Cloud can only exist in a narrow strip of the available pa- 
rameter space defined by the initial dark matter subhalo 
mass and the initial gas mass. This remains true with 
and without dark matter stripping. This narrow strip 
also appears, albeit shifted with respect to the present 
axes, if the Galactic halo is considered to be 10 x over or 
underdense from our model setup, or if the dark matter 
containing subhalo was virialized at an earlier redshift 
and correspondingly is denser than a subhalo virialized 
today. 

Consistent with our picture, the observed wake in Fig. 
1 presumably resulted from ablation processes due to the 
impulsive shock of the disk transit. We note with inter- 
est that the wake flares (i.e. expands transverse to the 
direction of the wake) beyond the putative dark matter 
core (Fig. 1). In our picture, this can be understood 
in the context of a transsonic flow. The Smith Cloud is 
moving at up to 300 km s -1 with respect to a hot halo 
with adiabatic sound speed ^200 km s _1 making the 
cloud-halo interface mildly supersonic. Once the stripped 
gas is left behind, it must undergo mixing with the hot 
halo through Kclvin-Hclmholtz instabilities, causing the 
stream to expand perpendicular to the flow. The temper- 
ature of the HI gas is predicted to rise along the wake. 
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Fig. 5. — The initial gas mass vs. the surviving tidal mass for three different dark matter models. The left hand column shows models 
which will have already been reduced to the final tidal radius by the beginning of the simulation; the right hand column shows models 
where tidal stripping is occurring at the same time as ram pressure stripping, simulating a cloud infalling for the first time. The allowed 
solutions for the Smith Cloud are shown in white. The upper and lower limits are discussed in j(5] 



Unless the confining core-dominated, dark matter halo 
is very massive, we do not expect the declining poten- 
tial seen by the wake material to influence the degree of 
flaring. These issues are being addressed in new hydro- 
dynamical simulations of the Smith Cloud that we will 
present elsewhere. 

With respect to further constraints on shock interac- 
tions, we expect that shock signatures at UV to x-ray 
wavelengths will have largely faded away, and the HI 
"hole" at the crossing point will have been substantially 
stretched by differential shear. For cloudlets smaller 



than 100 pc, thermal c onduction due to the halo corona 
(|McKee &: Cowiel 119771 ) and the Galactic radiation field 
convert the ablated gas to a clumpy plasma. We there- 
fore strongly encourage more extensive mapping of the 

HI and the ionized gas. 

As we discuss elsewhere (|Bland-Ha wthorr] [2008) . it is 
very difficult to see how the Smith Cloud, like several 
other large HVCs, could have come in from, say, the dis- 
tance of the Magellanic Stream. This was our motivation 
for considering a cloud stabilised against disruption by 
a confining dark matter halo. In essence, we have con- 
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Fig. 6. — The change in tidal dark matter and baryonic mass 
for a model over the last 150 Myr with input parameters (M v i r = 



1.2 X 10 9 M [M tid 



3 x 10 s Mq], M gas =l.l x 10 7 M Q ). The 



dark matter subhalos begin the simulation with a tidal mass of 
approximately 3 X 10 s Mq and a hydrogen mass of 1.1 X 10 7 Mq. 



sidered the prospect that the Smith Cloud constitutes 
a 'dark galaxy' where star formation never took place. 
Any evidence of a co-moving stellar population would 
have profound consequences. An essential requirement 
for star formation to occur appears to be the pr esence 
of molecular gas (e.g. iCalzetti fc KennicuttJ 12009). Wc 
envisage much deeper searches for the presence of molec- 
ular gas in compact HVCs, in particular, with the newly 
commissioned Cosmic Origins Spectrograph on the Hub- 
ble Space Telescope. 

Our model is not unique. The cloud may have been dis- 
lodged from the outer disk or confining potential of an 
infalling dwarf galaxy. A cloud metallicity of [Fe/Hjw-1 
may be appropriate in either scenario but there are no 
obvious candidates at the present time. The interloper 
must be on a prograde or bit which rules out som e in- 
falling dwarfs (e.g. wCen; [Bekki fc Freeman! 12003D but 
conceivably implicates disrupting dwarfs like Canis Ma- 
jor, assuming these were still losing gas in the recent 
past. 
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Fig. 7. — The change in baryonic to tidal dark matter mass ratio 
for a model over the last 150 Myr with input parameters (M v i r = 
1.2 x 10 9 M [Mtidal ~ 3 x 10 s M©], Mg as =l.l x 10 7 M Q ). The 
dark matter subhalos begin the simulation with a tidal mass of 
approximately 3 X 10 s Mq and a hydrogen mass of 1.1 x 10 7 Mq. 
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